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Abstract — We present a novel device-free stationary person 
detection and ranging method, that is applicable to ultra-wide 
bandwidth (UWB) networks. The method utilizes a fixed UWB 
infrastructure and does not require a training database of tem- 
plate waveforms. Instead, the method capitalizes on the fact that 
a human presence induces small low-frequency variations that 
stand out against the background signal, which is mainly affected 
by wideband noise. We analyze the detection probability, and 
validate our findings with numerical simulations and experiments 
with off-the-shelf UWB transceivers in an indoor environment. 

Index Terms — Device-free localization, UWB, indoor position- 
ing, signal processing. 

I. Introduction 

WIRELESS localization and tracking has attracted a 
great deal of research interest from the research com- 
munity, as location-awareness is fast becoming an essential 
feature in many application areas jl]. For indoor scenarios, 
ultra-wide bandwidth (UWB) transmission is a promising tech- 
nology, due to its high-resolution ranging and obstacle pene- 
tration capabilities |l2l-||6). Most practical UWB localization 
systems rely on targets (e.g., objects, people) to carry an active 
UWB device, which is used to facilitate time-difference-of- 
arrival or time-of-arrival measurements [0. In some scenarios 
(e.g., intruder detection, elderly care, smart environments, 
emergency response) it is desirable to have the ability to track 
people and assets in a passive manner, without requiring them 
to be equipped with any radio-frequency (RF) device. This 
is known as device-free localization, an overview of which 
can be found in the survey papers Q, 111]- Traditional device- 
free localization techniques were vision-based, relying on 
infrared motion detectors and video camera surveillance, but 
are limited to visible line-of-sight (LOS). Modem techniques 
overcome this problem through RF-based transmission, where 
received RF-signals are affected by the presence of people 
or assets in a quantifiable way. Research in this area can be 
broadly differentiated based on the narrow-band and wide- 
band nature of the signals involved. 

The earliest works in narrow-band device-free localiza- 
tion ll9l- lfm exploited human body-induced shadowing on 
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the transmission links by measuring received signal strength 
indicators (RSSI). In particular, (|9] detects motion and es- 
timates velocities based on variations in RSSI, while ifTOl . 
ifTTl employs filters or thresholds to detect changes in the 
RSSI. More recent works lfT2l . llT3l provided improvements, 
by partitioning the environment into cells and employing 
discriminant analysis in each cell 1121 and non-parametric 
statistical anomaly detection |fT3l . In lfT4ll - |fT8ll . received signal 
strength (RSS) variance is considered in an effort to improve 
robustness against environmental changes. Aspects of mobility 
on RSS variance were considered in llT6l - lfT8l . while models 
allowing distinction between mobile and static people were 
treated in |fT9l . 1201 . and for single and multiple people in 
121)1 . The term radio tomograpliic imaging was introduced 
in fT9l , as a general description for these systems, due to 
their similarity with computed tomography. While the low- 
cost implementation of narrow-band radios is quite attractive, 
susceptibility of the system to multi-path fading makes it hard 
to develop accurate models for dense, cluttered environments. 
These drawbacks can be overcome by considering larger RF 
bandwidths. 

UWB was demonstrated as an effective technique for 
human-being detection through respiratory movement in ll22l . 
l23l . Recently, experimental demonstrations were also given 
for different positions of the body and different antenna 
polarizations jm, MIMO UWB 123, and sensing of the 
person through obstructions ll26l . ||271 . For the estimation of 
respiration and heart rates, an analytical framework and a 
frequency domain technique were developed for single person 
in l28l . experimental results were given for multiple targets in 
l29ll . and related Cramer-Rao lower bounds were calculated in 
l30l . A new time-variant channel response model is introduced 
for breathing detection and human target ranging in fi7\, and 
the harmonics and the intermodulation between respiration and 
heart signals are analyzed in lISTl . While these works dealt 
with the detection of static people from breathing information, 
human-body detection and tracking were also studied exper- 
imentally for moving people in an open area in ll33l . 11341 . 
Passive object detection and tracking was also studied in UWB 
networks, in which an initial study on tracking was performed 
in l35l . which derived Cramer-Rao lower bounds, assuming 
a specular reflection model. Imaging of environments and 
objects based on a single UWB transmission was considered 
in l36l . l37l . and extended to multiple receivers in ||38l - PT]| . 
This has resulted in a flurry of research papers, investigating 
the impact of system geometry ||3FI , optimum detection met- 
rics l39l , imaging methods BTl . and signal analysis methods 
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P2l . While UWB device-free localization is able to overcome 
many of the drawbacks of narrowband signals, a number of 
issues are still unresolved. First of all, although experimental 
demonstrations were given for detecting the presence of people 
via UWB transmission, there is no unified signal processing 
technique for device-free detection, ranging, and localization 
of people. Second, in most existing methods, there is a need to 
develop a training database prior to localization. This database 
needs to be updated frequently, as changes in the environment 
may occur, precluding fast deployment. 

In this paper, we develop a novel device-free indoor human- 
body detection and ranging method, that can be applicable 
to UWB networks and used as a complementary method 
to existing device-free or device-based localization methods. 
Our detection method operates in the time domain unlike the 
works ll28l . 1291 . If32l and 1431 . in which frequency domain 
and polarization-based techniques are considered, respectively. 
Our method is based on the observation, corroborated by 
measurements, that a stationary human body introduces small 
temporal variations in the received UWB signal, even when 
the person is standing still. As such, it does not require any 
trainingQ Since the multi-path UWB signal, reflected from 
other objects is stationary over relatively long time windows, 
the deployment of this system in dense, cluttered environments 
does not affect the performance. Based on this observation, we 
describe a corresponding UWB signal and develop a generic 
detection technique. This detection technique is validated 
by numerical simulations and experimental measurements. 
Finally, we provide a limited set of localization results, based 
on experimental data. 

The remainder of this paper is organized as follows. In 
Section [III we formulate the problem. In Section |III] we 
describe a detection criterion, and analyze the corresponding 
false alarm and missed detection probabilities in Section |IV] 
Simulations and experimental results are provided in Section 
fVl before we draw conclusions in Section |VT] 

II. Problem Formulation 
A. Localization System 

We consider a system with N^. UWB radios (called beacon 
or anchor nodes), with a priori known positions e and 
a passive, human target with an unknown position x G K^. 
We assume that the anchors exchange signals which reflect 
off the target. Based on such a reflected signal from anchor i 
to anchor j, we obtain an estimate of the time of flight, or a 
corresponding distance 
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Fig. \. Mean of 100 UWB measurements, taken continuously within 20 
seconds in an indoor LOS environment, in the presence and absence of a 
person. 
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Fig. 2. Variation of the received waveform over time at delay instant 9.5 ns, 
in the presence and absence of a person. 



where the summation goes over all pairs of transmitters i 
and receivers j that have a measurement. Our goal is to 
determine di,j based on human body-induced signal variations, 
and without a template waveform. 



where e^j is measurement error. This enables the system to 
approximately locate the target on an ellipse whose foci are 
Xi and Xj, and with the length of the major axis equal to di,j. 
Given a collection of such measurements, the least-squares 
estimate of x is 

, 2 



arg min 



- X - X, 



(2) 



'in practice, some calibration may be required to detennine the ambient 
noise power and to set the system operating parameters appropriately. 



B. Signal Model 

In this section, we will develop a signal model that captures 
the effect of a human presence on UWB signals. Let us 
consider the following scenaricH: a transmitter and a receiver, 
separated by a distance of five meters in an indoor LOS 
environment, and a person on an ellipse around transmitter 
and receiver with major axis length 7.7 m (or, equivalently, a 
delay of 9 ns with respect to the direct signal path). Fig. [T| 

detailed description of the expeiimental setup is given in Section IV-BI 
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shows received UWB waveforms when the person is present 
(dashed Hne) or absent (solid line). We observe that the person 
affects the signal to some small extent around 9 ns after the 
first arrival. However, without a template (in this case the 
solid line) it is unclear how to determine the delay related 
to the person. A different view is offered in Fig. |2] showing 
the evolution of the signal at a fixed delay of 9.7 ns, in the 
presence and absence of a person. Such a signal is obtained 
by receiving a periodic waveform, and aligning the copies 
in the delay domain. As one would expect, in the absence 
of a person, the signal is essentially constant (here around - 
3500 analog-to-digital converter (ADC) counts) and affected 
by noise. Interestingly, in the presence of a person, the signal 
evolves slowly over time, due to minor temporal variations, 
e.g., induced by breathing. 

Based on these observations, we pose the following model. 
We transmit N^cp copies, with a period T^cp, of a ranging 
signal ID, 



sit) 



Ni-l 



(3) 



where Es denotes the energy of the signal, A^f is the number of 
pulses in s{t), bj G {±1} denotes the polarity code, Tf is the 
duration of the frame and p{t) denotes the unit energy UWB 
pulse of duration Tp < Tf, Cj € {0,1,..., Mi - 1} is the 
time-hopping code, Tc and N}^ are the chip duration and the 
number of chips per frame, respectively. The total duration of 
the signal is Tg = A^fTf, being smaller than the period T^cp. 
The receiver coherently combines A^f pulses during each of 
A^rep repetitions, leading to the following received signal 



''rop(^) = X] X! "'^(^ -Tl- fc^lcp) 



fe=l 1=1 

Nyap M 



^X Yap,i{t)p{t - Tpj - /cTrcp) + n{t), 



k=l 1=1 



(4) 



where the first term comprises L signal paths, being irrespec- 
tive of the presence of the person, in which ai and r; are 
the channel coefficient and the delay of the l-th discrete path, 
respectively. Assuming the delays are ordered, ti corresponds 
to the direct-LOS path. The second term in corresponds 
to the person, where in view of our model, there are M 
channel coefficients ap,i{t) varying slowly overtime due to the 
effects of the human body. Among them, ap^i{t) represents 
the signal path that is directly reflected off the person (i.e., 
direct reflection), whereas the remaining ones are related to 
signal paths that are reflected off the person and the other 
reflection sources in the environment (i.e, indirect reflections). 
Note that the delay associated with the direct reflected path 
Tp^i is always smaller than the delays associated with the 
indirect reflected paths Tp,; for I ~ 2, . . . , M, since the direct 
reflected path always travels a shorter distance. The noise 
n{t) is assumed to be zero-mean independent and identically 
distributed Gaussian noise with power spectral density No/ 2. 

Direct sampling of rrcp(i) at a sufficiently high rate W 



over a window Tdciay and aligning the N^^p copies, we 
obtain a two-dimensional model with r{kTi.cp,m/W) = 
rrcp{kT,cp + m/W), ap^i{kT,cp,m/W) = ap^i{kT,cp + 



m/W) and n{kT,^p,m/W) 
1,2, ... , Ni-cp and m = 0, 1, . 



= ri(fcr,.op + m/W) for k 

.,TdelayVF-l: 



r{kT,,p, m/W) 



L 

E 

1=1 

M 



aip{m/W " Ti) 



with 



/e; X apj(fcT,ep, m/W)p{m/W 
1=1 

n{kT,,p,m/W), 



E {n{kT,,p, m/W)nik%,p, m' /W)} 
NqW , 



-5{k~k')5{m-m'), 



(5) 



(6) 



where E{.} denotes the statistical expectation and S{.) is the 
discrete delta function. In view of our model, the second term 
in ^ again corresponds to the person, which introduces M 
channel coefficients that vary slowly over time (i.e., as a func- 
tion of k). We can further break up the signal corresponding 
to the person in a static part (i.e, independent of k) and a 
slowly-varying part as 

^/E',ap,i{kT,,p,m/W)p{m/W - Tpj) = 

^fE,C{m/W - Tp.O + ]{Tp,{)w{kT,,p)V(m/W - Tp.i) 

(7) 

where for a fixed m, C{m/W — Tpj) is a constant, depending 
on the propagation effects of the human body and the received 
signal power, /(rp.;) is a decreasing function in Tp_; to model 
the effect of propagation loss on the time-varying signal part, 
and 

/2(Tp,0 = E^GiTpj), (8) 

where G(rp_i) refers to the path gaiiH and defined as 



G(Tp,z) ^ G(r,ef) 



TpJ 



(9) 



for a reference path gain G(ricf ) at a reference delay r^cf and 
a path-gain exponent rj. The random variable w{kTi.cp) is of 
low frequency, zero mean and unit energy in time dimension 
(i.e., E {J2k ^'^^(^'^rop)} = -^rcp) and has unknown statistical 
properties in the delay dimension. Furthermore, w{kT,:cp) is 
assumed to be independent from the noise n{kTrcp,m/W). 
Finally, V{m/W ~ Tp,;) represents the influence from the 
transmitted pulse p{m/W — Tp,;) to the time-varying received 
signal part. 

Note that in our setting, L, a;, and t; are not known to the 
localization system, but assumed to be constant as a function 

'The distance dependence of the received signal power is typically ex- 
pressed by path-gain models. Given a reference received signal power P^.q at 
a reference distance dr.o and path-gain exponent r), the received signal power 
at distance d is usually expressed in dB as Pr = P^ q + 10»7logjQ(c(r.o/d). 
Based on our experiments, we have chosen a similar model to show the 
relation between the energy of the slowly-varying signal and the the delay of 
the reflected paths. 
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of time, during 1 < A: < A'l-cp- Moreover, for the transmission 
between transmitter i to receiver j, the human-body-induced 
delay of our interest is given as 



Tl 



(10) 



where c is the speed of light. Since we can estimate ti using 
standard time-of-arrival estimation techniques [[3], and since 
we know the positions of the anchors, and thus also | |xi — Xj 1 1, 
from ( [Tol l it is straightforward to determine di,j, defined 
in (HJ. Furthermore, we will abbreviate r(fcTrcp, m/VF) as 
r(fc, m) with the understanding that samples in the different 
dimensions (time, delay) are taken at different rates. In this 
proof-of-concept study, we will not consider the statistics of 
the measurement error e^j, or the impact of mobility and 
corresponding tracking methods. 

Our main objectives are to determine (i) if the person is 
present and (ii) on which ellipse the person is (i.e., what is 

^P.i?)- 

Comment: We should point out that our signal model as- 
sumes no other sources can create low-frequency disturbances 
to the UWB waveforms. In practice, slowly moving equipment 
or periodically moving stationary objects (e.g., a fan) may have 
similar effects as a human being on the UWB waveforms. 
Further study is required to incorporate these error sources. 

III. Detection and Device-Free Ranging 

In this section, we develop a device-free person detection 
and a ranging technique, based on the signal model in (|5]l. We 
will proceed as follows: we first consider a fixed delay m/W 
and determine a meaningful (though not necessarily sufficient 
or minimal) statistic. Then, we combine these statistics over 
multiple delays, allowing us to detect the presence of a person 
and to infer Tp i. An analysis of the final statistic is deferred 
to Section HV] 

A. Statistic for Fixed Delay 

Considering a fixed delay m/W in (|5]l, we obtain signals 
as shown in Fig.|2] As the mean of the signal does not convey 
relevant information for our purpose, we will subtract it, thus 
leading to the following signal model 



rimik), no person affects delay m 

Xm(k) +n,nik), person affects delay m, 

(11) 

where we have moved the delay index m to a subscript to 
emphasize the dependence on the time dimension k. Here, 
Xm{k) is a low-frequency signal induced by the presence of 
the person and nm{k) is assumed to be white Gaussian noise 
with variance NqW/2. In order to filter out high-frequency 
noise components from the low-frequency variations induced 
by the person, we decompose the received signal as 



rmik) = rrn,h{k) + r„i,H(fc) 



(12) 



where r„i^i^{k) and r,„_H(fc) are the low- and the high- 
frequency components, with fractional bandwidths (i.e., the 
bandwidth that is normalized to l/Trop) of /3 and 1 — /3, 
respectively, where f3 <C 1. In practice, f3 depends upon 



how fast the movements of the person occur within the time 
window of TrcpA'rcp- For instance, fS/T^cp should be around 
0.2-0.4 Hz for normal breathing. Stacking the time samples 
yields the vectors r„i, x™, and r„i l- A likelihood ratio test 
then yields a test statistic 



exp I 



NoW 



2r^ 



exp 



NoW 



exp 



NoW NoW 



(13) 



Since x^ is unknown, we treat it as a deterministic vector 
in the likelihood ratio test and replace it with an estimate. 
We only exploit the low-pass nature of x^, and thus set the 
estimate to x,„ = r„i i^. We easily find a final statistic, after 
normalization with /3A^rop afl 



y{m) = 



\r,n.U 



(14) 



Nol3WN,,p 

where, in the absence of the person, E{?/(m)} ~ 1. 

B. Statistic for Delay Window 

Up to now, we focused on a single delay value, however 
the person will affect a window of delays. Assuming that the 
person has an effect over the duration of the transmitted pulse 
Tp, we can aggregate the information over multiple delays by 
averaging the delay-specific statistic over a window around a 
trial delay r, which is assumed to be an integer multiple of 
1/W: 

iT+Tj2)W 



D[r) 



1 



T^W 



E 

m=(T-Tp/2)VK 



y[m). 



(15) 



The presence of a person can thus be determined by comparing 
D{t) to a threshold. Hence, the person is detected as 



D{t) < 7 no person present 
D{t) > 7 person present, 



(16) 



where 7 is a threshold. The selection of 7 depends on the de- 
sired performance trade-off, and will be treated in Section |IV] 

In addition to the presence of the person, D{t) also conveys 
information about the delay of the human body induced 
reflections, from which we can obtain an estimate of Tp 1. 
Here, we introduce three practical estimation criteria: 

1) Line Search: This criterion is based on the selection of 
the largest sample in D{t), and specifically given by 



To 1 = are max D(t) 

re[0,Tdelay] 



(17) 



2) Threshold Crossing: The maximum value of D{t) needs 
not always correspond to Tp 1 but may come from a later 
arriving indirect reflected path Tp./ with I ~ 2, . . . , M. A 
detailed discussion of this phenomenon is left to Section [V-BI 
To provide robustness against these failures, we introduce 
another criterion where we estimate Tp.i based on the first 

"^In practice, A^o must be estimated as well. This can easily be achieved 
using standard techniques 1441 . We will return to this in Section FV-BI 
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threshold crossing: 



Tp,i = mill {r : D{t) > 7} 



(18) 



7 is also a threshold and can be different than 7, depending 
on the desired performance and false alarm criteria. 

3) Maximum Rise Search: Although the previous approach 
provides additional robustness, its performance is dependent 
on the threshold value, which may be hard to determine in 
some cases. We can also exploit the shape of D{t) in which 
we search for the maximum increase and thereby remove the 
requirement for the threshold determination. Specifically, the 
maximum rise search criterion is given by 

^(T + Twi„/2)-Z?(r-rwin/2) 



Tp 1 = arg max 

re[0,Tdoi 



T 

J- W 



(19) 

where Twin is an arbitrary window (whose length we will set 
to Twin = Tp, neglecting the pulse distortion in the channel). 

IV. Performance Analysis 

In this section we determine the false alarm and missed 
detection probabilities of the detector proposed in (IT&t , for an 
estimate r*, which can be generated from any of the three 
above-mentioned criteria. 

A. Probability of False Alarm 

In the absence of the person, substitution of (fl4l i into ( fTSl l, 
and accounting for the fact that rm,L(fc) = nm,hik) is low- 
frequency noise, yields 

iT'+Tp/2)W 

^(-^^ = TpNom^N..p ^ E (20) 



For sufficient number of time samples A^tcp, we can invoke 
the central limit theorem (CLT), and approximate D{t*) ^ 
M{fJ.f, erf), where 



fi{ = 1 

and CTf is derived in Appendix |A] and found to be 



' TpPWN,,p 
Hence, the probability of false alarm TpA is given as 

TpA = Pt{D{t*) > 7 1 110 person present} 
'7_/if ' 

(Tf 

where Q{.) denotes the Q function. 



(21) 



(22) 



(23) 



B. Probability of Missed Detection 

In the presence of the person, substitution of (fl4] i into (fTSl l, 
and accounting for the fact that r„i^L(fc) — Xm{k) +nmx{k), 
yields 

(r*+Tp/2)iy 

Dir*) 



1 



TpiVo/3W^2^rcp 



E 



m={T*-Tp/2)W 



(24) 



From the right-hand side of Q, is written for the direct 

reflected path as 

Xraik) = fiTp^l)w{k)V{m/W - Tp^). (25) 

Substitution of (l25T l into (l24b . expanding the square, and 
considering that V^{m,/W — Tp^) — W yields 



D{t*) = 



TpNom^N,,p 



{T'+Tp/2)W 

E ( 

m = {T*-Tp/2)W 



\xm\\^ + ||n„, lII^ 



G(rp,i)g, 
TpNol3W^N,,p 



TpNopW^N, 



rep 



TpNopW^N,,p 



(T*+Tp/2)W 

E V^m/W - rp,i) \\^ 

m={T*-Tp/2}W 
(r*+Tp/2)H' 

m={T*-Tp/2)W 
(r*+Tp/2)W 

x^^n„,L- (26) 

m={T'-Tp/2)W 



Invoking again the CLT for sufficiently large A'rep, D{t*) 
JV{iJ.d,crl), where 



^ G(Tp^ 

TpNoPW 
and (see Appendix |B] for details) 



(27) 



2 



+ 



AGiTp,i)E, 



TpPWN,,p T^Nol3-^Wm,,p 

Finally, the probability of missed detection Pmd is obtained 

as 



^'md = Pr {D{t*) < 7 I person present} 

'7 - Md' 



1-Q 



CTd 



(28) 



V. Performance Evaluation and Discussion 

In this section, we quantify the performance of the detection 
and the device-free ranging technique from Section |lll] We 
will first consider the detection and provide numerical results 
for the relevant performance measures. Then, we will show 
experimental results for both detection and the device-free 
ranging. Finally, the localization performance will be evalu- 
ated. 

A. Detection and Device-Free Ranging: Simulation Results 

Simulation Setup: Unless otherwise stated, we consider the 
7-th derivative Gaussian pulse as a transmitted signal with a 
duration of Tp = 1.4 ns and set the receiver sampling rate to 
61 ps, resulting in 23 delay samples over the transmitted pulse 
duration. We collect a varying number of repetitions with a 
period of Tiop = 0.2 s. Without loss of generality, the variation 
due to the person is considered to have a sinusoidal shape ll28l . 
ll30l with a normalized frequency of 0.04 cycles/sample (i.e., 
an absolute frequency of 0.2 Hz). The fractional bandwidth 
of the filter, applied over time, is set to /3 = 0.1 (i.e., an 
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Fig. 3. False alarm probability versus threshold under different number of 
samples over time. 
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Fig. 4. Missed detection probability versus threshold under different number 
of samples over time at low (15 dB) and high (30 dB) SNR and do = 6 m. 



absolute bandwidth of 0.5 Hz). Hence, it is fully possible to 
capture the temporal variations induced by the person. Based 
on our measurement datc|f|, we set the model parameters in 
^ as G(rrof) = 3.6 x IQ-^ r,.of = 17.8 ns and ^ = 5.5. 
Finally, the device-free range is denoted as do — cTp i, and 
the received signal-to-noise ratio (SNR) is defined as SNR — 
Es/Ni^, where we have scaled the signal with 1/ai. 

'We first define a reference received signal as the first TpW samples of 
the received wavefonn after the leading edge at a given transmitter and the 
receiver distance. For the same distance, we calculate the total energy of the 
human-body induced delay samples (again over the length of TpW) after 
removing the mean from the signal. Finally, we divide by the energy of the 
reference signal. While extracting these parameter values, we only consider 
the cases when the person is not blocking the direct line-of-sight (i.e we 
exclude the shadowing effect induced by the person). 
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Fig. 5. Missed detection probability versus SNR under different false alarm 
probabilities (do = 6 m, A^rep = 1000). 



Results and Discussion: We first validate false alarm and 
missed detection probabilities derived in the previous section. 
The validation was based on simulating the system described 
above and determining the false alarm and missed detection 
probabilities through Monte Carlo estimation. Fig. |3] shows 
the false alarm probability as a function of the threshold 7 
for different values of A^icp- Note that this performance does 
not depend on the SNR. We observe good agreement between 
the simulations and the predicted performance, especially for 
larger values of iViop, as involving the CLT becomes more 
valid when more observations are collected. Fig. |4] shows the 
missed detection probability as a function of the threshold 7, 
for two different SNR values, and different values of A^tcp- 
We again observe good agreement for large iVrep, for both low 
and high SNR. When A^icp is low, the predicted performance 
is worse than the simulated performance. This is due to the 
fact that the Gaussian approximation does not completely fit 
to the distribution of ( l24l i. The mismatch occurs in left-side 
tail, where the Gaussian approximation underestimates the true 
probability density. Note that, in general, A'jop is dependent 
on both total observation duration over time window and the 
measurement update rate (l/Trcp). Therefore, the performance 
can be described well with analytical results for high update 
rates over shorter time windows (sufficiently large to capture 
slow variations) or slower update rates (again sufficiently fast 
to capture slow variations) over longer time windows. In the 
following part, we will only show the analytical results assum- 
ing that update rate is sufficiently high to collect A^icp — 1000 
samples over the time window. 

Figs. |5] and Fig. |6] show the dependence of the missed 
detection probability on SNR and the distance, respectively, 
for varying values of the false alarm probability PpA- Fig. |5] 
illustrates that for a given device-free range do and Pfa, we 
need about 1 dB of SNR improvement for every order of 
magnitude reduction in Pmd (e.g., around 3 dB to get from 
Pmd = 10^^ to Pmd = 10~^). Furthermore, the performance 





Fig. 8. UWB setup for device-free localization measurements, the person is 
standing at HI with a perpendicular body orientation. 



Fig. 6. Missed detection probability versus device-free range under different 
channel conditions (PpA = 10"^, A^rep = 1000). 
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Fig. 7. The overview of measurement scenarios: the transmitter (Tx), the 
receiver (Rx) (red cross) and the subject positions (blue rectangular) and 
device-free detection range (dotted ellipse) limited by the current radios. 



of the system also depends on the position of the person 
with respect to the transmitter and the receiver distance, as 
we observed, the slow-variation of the signal attenuates by 
the increase in the reflected-path length. Fig. |6] shows that, 
considering Pmd = 10"'^, it is possible to detect the person 
up to 10 m device-free range at high SNR (e.g., 35 dB). Note 
that the device-free range is always larger than the arrival 
time of the signal (do > ^i)- Therefore, although results show 
that at higher SNRs the device-free range increases, the SNR 
decreases with the distance between the transmitter and the 
receiver due to the path loss. Thus, in practice the device-free 
range is limited by the transmitter and the receiver distance. 



B. Detection and Device-Free Ranging: Experimental Results 

Experiment Setup: We performed experiments in a fitness 
room at the Chalmers University of Technology building 
with two identical, commercially available and FCC-compliant 
UWB radios (Time Domain P400 radio^. The radios are ca- 
pable of performing communications and ranging using UWB 
signaling within the frequency range of 3-5.5 GHz. They were 
equipped with identical antennas that are omni-directional 
within the operating frequency range. While these radios are 
primarily designed for two-way time-of-arrival ranging, they 
are capable of capturing a waveform. At the receiving radio, 
the waveforms are sampled with a sampling period of 61 ps 
over a window of approximately 10 ns, limiting the maximum 
captured multi-path length to be less than approximately 3 m. 
With the current hardware release, beyond the 10 ns window, 
captured waveforms are no longer stable and processing across 
time is not meaningful. The transmitted pulse duration was 
measured with high-frequency oscilloscope and found to be 
approximately 1.4 ns. 

Measurements were performed with fixed transmitter (Tx) 
and receiver (Rx) positions, while the person was standing at 
13 different points as shown in Fig |7] For each position, the 
person stood in two orthogonal orientations, so that the body 
is either perpendicular or parallel to the line-of-sight between 
radios. During the experiments, the radios were mounted on 
tripods at a height of 1.23 m above the ground and were 
connected via ethernet cables. Measurements were taken using 
the graphical user interface at a rate of 5 measurements/second. 
The waveforms were collected over 20 seconds, allowing us to 
obtain 100 snapshots over time for each scenario. During the 
experiments, care was taken to keep the environment static. 
Therefore only the subject was standing within the device- 
free detection range limited by the radios. The operator, who 
was controlling the experiments, and all the other people 
were outside the range. A typical measurement scenario is 
shown in Fig. [8] During the off-line postprocessing, fine 
alignment around the leading edge point was performed by 

^See |http ://www.timedomain.com7| 
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cross-correlation to make sure that all the waveforms are 
aligned over the time window. The starting instant of the 
waveform (the zero delay instant) is defined based on the 
leading edge detection point, provided by the radio. Similar 
to the simulations settings, we applied the low-pass filter with 
a fractional bandwidth of (3 ~ 0.1. We chose the detection 
threshold as 7 = 1.38 to obtain a false alarm probability of 
PpA ~ 10~^ (see also Fig. |3]l. 

Note on background noise and tuning jitter: The radios are 
subject to a number of hardware impairments, in particular 
the timing jitter In general, the timing jitter causes deviation 
of the transmitted pulses from reception at integer multiples 
of sampling time ||45l . At a certain delay sample, the signal 
amplitude varies over time proportional to the time rate 
of change of the received waveform (i.e., the derivative of 
received waveform to time). This variation wiU be greater for 
higher amplitudes, since the change in the amplitude will be 
greater, leading to a higher impact of the timing jitter for these 
samples. This resulted in different amounts of background 
noise for delay with high signal values, compared to delays 
with low signal values. For that reason, we estimate Nq on a 
delay by delay basis. For a given delay, m, we estimate the 
noise power No^m as 

where rm,H(fc) = rm{k) — rm.hik), is the part of the signal 
that contains only noise, irrespective of the presence of a 
person. The estimate was then substituted in (fT4] i to allow 
detection and ranging in the presence of hardware imperfec- 
tions. 

Experiment Results and Discussion: Experimental results 
for the device-free ranging are tabulated in Table |l] showing, 
for each of the measurement locations, the true device-free 
range, the device-free range error (calculated as the difference 
between the estimated and the true device-free ranges, and 
averaged over the body orientations) for three ranging criteria 
introduced in Section lTlI-BI Our results show that the proposed 
method was able to detect the presence of the person in all 
cases. Furthermore, when we employ the line search approach, 
the device-free range error is found to be always positive 
including different orientations and less than 60 cm for all the 
cases, except the case where the person blocks the direct line- 
of-sight (case H4). A deeper investigation revealed that tempo- 
ral variations are not only observed at a single delay window, 
but at multiple delay windows of the received waveform. This 
effect is not dominant in other cases and may be due to the 
fact that the signal is reflected by another reflection source, 
after being reflected by the person. We also observe that the 
estimated noise power may be lower for samples related to 
these indirect reflections, resulting in higher values for the 
decision statistic (this is a side-effect of the way we dealt with 
timing jitter). On the other hand, when we employ threshold 
crossing approach, we obtain a better ranging performance as 
shown in Table |T] for position H4. In this case, we kept the 
threshold 7 = 7 = 1.38. Although the ranging performance 
improves for H4, the threshold crossing criterion results in 



TABLE I 

Experimental Results for Device-Free Ranging 









Error [m] 




Position 


Range [m] 


Line 


Threshold 


Maximum 






search 


crossing 


rise search 


HI 


7.85 


0.2 


-1.22 


0.01 


H2 


6.43 


0.4 


-0.55 


0.09 


H3 


5.4 


0.55 


0.15 


0.17 


H4 


5 


1.4 


0.42 


0.81 


H5 


5.4 


0.6 


0.25 


0.29 


H6 


6.43 


0.38 


0.01 


0.05 


H7 


7.85 


0.19 


-1.29 


-0.12 


H8 


6.7 


0.27 


-0.85 


-0.14 


H9 


7.39 


0.27 


-1.07 


-0.03 


HIO 


7.5 


0.07 


-1.24 


-0.31 


Hll 


6.7 


0.44 


0.06 


0.02 


H12 


7.39 


0.32 


-0.13 


-0.04 


H13 


7.5 


0.46 


-2.07 


0.17 



early detections in some other positions. This is again mainly 
due to the variation of the estimated noise power across delay 
samples, which has lower values for early delay samples, 
resulting in higher values for the delay statistic. Increasing 
the threshold may lead to better results, though, in general, 
the determination of the optimum value might be a difficult 
task. Finally, as the maximum rise search criterion removes 
the requirement of the threshold determination, it outperforms 
the threshold crossing and gives a better ranging performance 
compared to the line search for H4. 

Note that in our current work, we did not investigate the 
effect of other moving objects on the detection and the ranging 
performance. The detection method might give false alarms 
because of the changes in signal due to other autonomous 
systems, hence reducing the applicability of the system in 
some environments. Such changes in signal might depend on, 
for instance, the main material of the object and its shape 
and/or the strength or character of the movement, and should 
be considered by inspecting the dynamic nature of the environ- 
ment when it is not populated (e.g., an industrial environment 
can be highly dynamic due to the machinery, whereas office 
environments are, most of the time, less dynamic). Effects due 
to the other non-stationary objects are left for future work. 

C. Device-Free Localization: Experiment Results 

While our current analysis did not include a statistical evalu- 
ation of the device-free range error, we here present indicative 
results of the localization capabilities of the proposed system. 
Experiments were performed in a room, located on the top 
floor of the Carre building in University of Twente. The room 
was partly furnished with tables and chairs, and also has thick 
metallic pipes for ventilation and water supply. We deployed 
four anchor nodes, connected to a computer via network 
cables and a switch, and arranged in a square at positions 
XI = (-5,0), X2 = (0,0), X3 = (-5,-5), X4 = (0,-5). 
A person stood on 24 different positions in a grid within 
a 3 meter by 5 meter area. For each position, we took six 
anchor-to-anchor measurements. Between each anchor pair. 
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Fig. 9. Experimental results for localization performance, obtained when the 
person stands on 24 different positions in an indoor environment with four 
anchor nodes (marked as red squares). Threshold crossing criterion is applied 
to obtain the device-free range estimates. 



measurements were performed over 20 seconds with a rate 
of 50 measurements/second (i.e., 1000 snapshots). 

In Fig. |9l locaHzation resuhs are shown when we apply the 
threshold crossing criterion for device-free range estimates. 
The threshold is chosen as 7 = 1.1, which corresponds 
to PpA of 10~^ for iVrcp = 1000 based on the numerical 
results from Fig. |3] We get error values ranging from 0.12 m 
(measurement point 18) to 1.8 m (measurement point 3). 
The root mean square (RMS) localization error, computed 



where A^pos denotes the 



as Y '^/NposJ2i=r ll^p,' ~ ^ 
number of standing positions of the person (i.e., 24 in this 
case) is 1 m. When employing the maximum rise search and 
line search methods, the RMS localization error increases to 
1.33 m and 1.6 m, respectively. Higher RMS localization 
errors can be explained similarly to our observations from 
Section IV-BI In general, significant negative errors in device- 
free range estimates with threshold crossing criterion (i.e., 
when the person is close to the direct-LOS between devices), 
and significant positive errors with line search and maximum 
rise search criteria (i.e., when the person is away from the 
direct-LOS), leading to increased errors in the final position 
estimates. 

VI. Conclusion AND Future Work 

In this paper, we developed a novel indoor UWB device- 
free person detection and ranging technique that does not 
require any knowledge about the environment. Our method 
relies solely on exploiting the temporal variations in the 
received signal induced by the presence of the person. We 
observed that the signal evolves slowly over time, even when 
the person is standing still. Based on this, we developed a 
signal model. Using this model, we formed a decision statistic 
to detect the person and estimate the delay of the human- 
body reflected signal from which we obtained the device-free 
ranges based on the three different criteria. We analyzed the 



detection probability and validate with simulation results. We 
also demonstrated the performance with experiments, which 
revealed that our method was able to detect the person and (for 
most of the cases) resulted in satisfactory device-free ranging 
performance. However, outliers in range estimates degrade the 
localization performance of the system which we also showed 
in an experimental environment. Future avenues of research 
include the development of ranging error models, locating 
multiple people, tracking of moving people and investigating 
the effect of other non-stationary objects in the environment. 
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Appendix 

A. Derivation of Statistics for False Alarm Probability 
We recall that in the absence of a person 

(r*+Tp/2)M/ 



D{t*) = 



1 



■ E 

m={T'-T^/2)W 



n„, l| 



(30) 



where nm.-L{k) are Gaussian noise samples obtained after low- 
pass filtering and with variance NqW (5. We first calculate 

1 



.{D\r*)] 



(TpA^o/3M^2^rcp) 
1 



5]E{||n„ 



ijj{m^ m!) 



(31) 



where the summation goes over interval of length TpW (in 
two dimensions) and where we have introduced 

1 



?/'(m, m!) = 



^E{<L(A:)n2„,_L(fc')} 



(32) 



When m ^ m', we find that 

ih{m,m') ~ — ^ — 7^ — t: — T7— 



(33) 



■°p k.k' 

= 1. 

While for m ~ m', we can make use of the fact that 

E {n^„x(fc)"™x(fc')} - E {<L{k)}E {n^x(fc')} 

+ 2(E{n™,L(fc)«m,L(fc')})' 
= {NoWpf + 2{NoWg{k~k')f, 

(34) 

where g{k) is a discrete low-pass filter with gain 1 in the band 
[~l3W/2, 13W/2] and where we have utilized the fact that 



K{n,nMk)n,nMk')} - N„Wg{k - k'). 



(35) 
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Hence, 

+ 2 {NoWg{k - fc'))'] 

\{NoWPN,,,,) 



We easily find that 



E{i?(r*)} = ^i^i^ + 1. 



(40) 



k,k' 



^ [(iVoM^/37V,ep)' 



CXD 

k— — oo 

1 



Recalling that odd-order moments of zero-mean Gaussian ran- 
dom variables are zero, we determine E{Z?^(r*)} explicitly 
as 

^ ^ \TpNol3W J 

E< I .L,. EE<L(fc)) 

ni k / 



y^y^nm,L(fc)a:m(fc) 



= 1 + 2/(/3iV^.ep) 

Substitution leads to 



(36) 



m k 



(41) 



i^(Tpl^x (l + 2/(/?7V,ep)) 



y 2 14/2 



Here the second term simplifies to 1 + 2 j (T^fiW N -rc^) , the 
fourth term to 2G{t^,x)EsI{T^'NqP,W). The third term is 
equal to 

Terms = 



Tp7Vo/3W^27v, 



rep 



E E'^{"'".L(fc)n™',L(fc')} 



m,m' k.k' 



(T^W^ - TpW) X 1 



"2^/2 
1 



xE{x^{k)Xra'{k')} 



2TpW/{m,p) + T^W^) 



TpNopW^N, 



rep 



= 1 



So that we finally obtain 



a! 



TpPWN,, 



(37) 



(38) 



'^A{m)+ ^ B{m,m') 

I m.m'^7n 

(42) 



B. Derivation of Statistics for Missed Detection Probability 
We recall that in the presence of a person, 

(r*+Tp/2)lV 

Dir*) 



in which 

A{m) = E {?i„t,L {k)n,n,h {k')}E {.t,„. (fc)a;,„. (fc^) } 

= iVoM^EEff^^ - fc')E{a;™(fc)x™(fc')} 
fc fc' 

= NoWY,^{xlik)}, (43) 



1 



rp7Vo/?Vl^2^rcp 



E (II-" 



l' + l|n™.Lf 



and B{m,m') — 0. Hence, 



Term':! = 



m=(T*-Tp/2)iy 



_ G(Tp,i)£;, 



(T*+Tp/2)Vy 

X E ^^(WW^-T-p,i) ||wf 

m=(r*-Tp/2)Vy 

(r-+Tp/2)W 
E 

m=(r*-Tp/2)W 
(r-+Tp/2)lV 



4jVoiyG(Tp,i)^, 

p OA" rop „j ^ 

4G'(rp,i)^, 



EEE{^™(fcf} 

(44) 
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jyVo^W^ep 



Putting everything together, we find 



Tp/3M^iVrcp 



Tp/^VFiVo V Tp/JW^A^r, 



(45) 



+ 



E """"'^ 



(39) 



m={T*-T^/2)W 
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Subtracting 

EMi^(r*)}^l + ^^i^+r^i^V, (46) 



gives us 



(47) 



2 / 4G(Tp^ 
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